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The fully nonlinear dynamics of spin and charge in spin-Calogero model is studied. The latter is an 
integrable one-dimensional model of quantum spin- 1/2 particles interacting through inverse-square 
interaction and exchange. Classical hydrodynamic equations of motion are written for this model 
in the regime where gradient corrections to the exact hydrodynamic formulation of the theory may 
be neglected. In this approximation variables separate in terms of dressed Fermi momenta of the 
model. Hydrodynamic equations reduce to a set of decoupled Riemann-Hopf (or inviscid Burgers') 
equations for the dressed Fermi momenta. We study the dynamics of some non-equilibrium spin- 
charge configurations for times smaller than the time-scale of the gradient catastrophe. We find 
an interesting interplay between spin and charge degrees of freedom. In the limit of large coupling 
constant the hydrodynamics reduces to the spin hydrodynamics of the Haldane-Shastry model. 
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I. INTRODUCTION 



One-dimensional models of many body systems have been a subject of intensive research since the seventies. Due 
to the low dimensionality, standard perturbative approaches developed in many body theory are often inapplicable. 
On the other hand some techniques specific to one spacial dimension are available and allow to treat systems of 
interacting particles non-perturbatively. The Fermi Liquid paradigm is replaced by the Luttinger Liquid theorjffl 
in one dimension. One of its most striking predictions is that at low energies spin and charge degrees of freedom 
decouple. One can say that at low energies physical electrons exist as separate spin and charge excitations. At higher 
energies it is expected that spin and charge recombine into the original electrons. One can see the traces of spin-charge 
interaction taking into account corrections to the Luttinger liquid model arising from the finite curvature of band 
dispersion at Fermi energjff'. The coupling between spin and charge in one-dimensional systems was studied both 
perturbatively and using integrable models available in one dimension^. 

In this paper we study the interaction betwee n spi n and charge in another integrable model ~ the spin-Calogero 
model (sCM). This model is a spin generalizatiorP'-'^ of the well-known Calogero-Sutherland model^. 

Calogero-Sutherland type models occupy a special place in ID quantum physics. They are exactly solvable (inte- 
grable) but are very special even in the family of integrable mo dels. In par ticular, they can be interpreted as systems 
of "non-interacting" particles with fractional exclusion statisticgSEEMloill], 

The sCM model is given by the following Hamiltonian: 

where we took the mass of particles as a unity and Pji is the operator that exchanges the positions of particles j and 
PI The ± sign in the exchange term corresponds to ferromagnetic and anti-ferromagnetic ground state respectively if 
we are studying fermions. Similarly, it corresponds to anti- ferromagnetic and ferromagnetic ground state respectively 
if we are considering bosonic particles. The four cases can be summarized as: 



Bosons 
Fermions 



Anti-ferromagnetic , 
Ferromagnetic , 

Ferromagnetic , 
Anti-ferromagnetic . 



The coupling parameter A is positive and N is the total number of particles. 

As it has been already noted above the sCM is a very special model. In particular, in contrast to more generic 
integrable or non-integrable models the spin and charge in sCM are not truly separated even at low energie^. 
Of course, one can still describe the low-energy excitation spectrum of sCM by two independent harmonic fiuid 
Hamiltonians, one for the charge and the other for spin. However, it turns out that for the sCM the spin and charge 
velocities are the sam^, i.e. spin and charge do not actually separate. 

Here we study the spin-Calogero model in the limit of an infinite number of particles using the hydrodynamic 
approach. Albeit, the colle ctive field th eory /quantum hydrodynamics of the spinless Calogero-Sutherland model has 
been studied in great detaiP^ESElIIIIill^ a complete understanding of its spin generalization is still lacking although a 
considerable progress has been done recently in Refs.'^^'^^. 

We study the nonlinear collective dynamics of sCM in the semiclassical approximation, additionally neglecting 
gradient corrections to the equations of motion. This limit is justified as long as we consider configurations with small 
gradients of density and velocity fields. The gradientless approximation is commonly employed in studying nonlinear 
equation^i^ and allows to study the evolution for a finite time while the first nonlinear contributions are dominant. 
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For longer times, the solution will inevitably evolve toward configurations with large field gradients (such as shock 
waves) and the gradientless approximation becomes inapplicable. Nevertheless, in the initial stage of the evolution, 
corrections due to gradient terms in the equations of motion can be neglected (see further discussion in the Sec. VB I. 
We derive the gradientless hydrodynamics Hamiltonian from the Bethe Ansatz solution of the model. 

The paper is organized as follows. In Sec. [ll] we start with the simplest spinful integrablc model - a system 
of free fermions with spin. We briefly review the Bethe Ansatz solution for spin-Calogero model in Sec. |III| and 
deduce the hydrodynamic Hamiltonian for the sCM from this solution in Sec. |IV| neglecting gradient corrections. The 
corresponding classical equations of motion are given in Sec. |V] It is shown that variables separate and the system 
of hydrodynamic equations is decoupled into four independent Riemann-Hopf equations for a given special linear 
combinations of density and velocity fields - the dressed Fermi momenta. In Sec. |VI| we illustrate that in the limit of 
strong coupling the hydrodynamics of sCM is reduced to the hydrodynamics of Haldane-Shastry lattice spin model 
giving the hydrodynamic formulation of the so-called freezing tricl¥^. We present some particular solutions of the 
hydrodynamic equations demonstrating nonlinear coupling between spin and charge degrees of freedom in the sCM 
in Sec. |VII| and conclude in Sec. |VIII| To avoid interruptions in the main part of the paper some important technical 
details are moved to the appendices and are organized as follows. In Appendix [X] we use asymptotic Bethe ansatz to 
derive the hydrodynamics of sCM and to explain why variables separate in this system. In Appendix [B] we describe the 
notion of true hydrodynamic velocities. In Appendix [C] we relate the hydrodynamics of sCM to two infinite families 
of mutually commuting conserved quantities and collect our results for the hydrodynamics in the different regimes 
of sCM. Finally, in Appendix [P] we derive a hydrodynamic description of the Haldane-Shastry model from its Bethe 
Ansatz solution. 



II. FREE FERMIONS WITH SPIN 



For one-dimensional free fermions without internal degrees of freedom the lowest state with a given total number of 
particles and total momentum corresponds to all single-particle plane wave states filled if the corresponding momentum 
k satisfies < k < fc^. Here kL^R are left and right Fermi momenta respectively which are defined by the given 
number of particles and momentum of the system: 



'^^ dk 
, 2^ 



27r 



— hk = h-^ 

2tt 4tt 



pv . 



(2) 
(3) 



N/L = 
P/L = /' 

Here we introduced the (overall) velocity of the system v which is given from (2j3l by 

,^ kR + kL , , 

v/h = . (4) 

Inverting (2]4| we express the left and right Fermi points k^^R in terms of the density p and velocity v as 

kR^L ^ v/h±Trp. (5) 

The energy of this state is given by 



E/L ■ 



dfc tfk^ ^ykl 



2n 



12n 



pv 



6 



(6) 



Up to this moment p, v, kR^^ £^re just numbers characterizing the chosen state of free fermions (only two of them 
are independent). Assuming the locality of the theory we promote these numbers to quantum fields and write the 
hydrodynamic Hamiltonian of free spinless fermions as 



n 



dx 



p{x)v'^{x) 



2^2 



~6 



■p\x) 



dx h 



2[kR{x)] - [kL{x)] 



127r 



(7) 



Here we consider p{x) and v{x) as quantum field operators of density and velocity (and kR^^ as given by ([5|)) having 
canonical commutation relations^l- 



[p{x),v{y)] = -iM'(a; - y) 



(8) 
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Of course, gradient corrections to ([t]) are generically present and the above "derivation" is just a heuristic argument 
(semiclassical in nature). It turns out that ([7| is, i n fact, exact for free fermionsl^ It can be derived rigorously either 
using the meth od of c ollective field theorji^^^El or conventional bosonization technique (but without linearization 
at Fermi points)l32ll3Zl. 

The two terms of (|7| have a very clear physical interpretation. The first term is the kinetic energy of a fluid moving 
as a whole - the only velocity term allowed by Galilean invariance. The second one is the kinetic energy of the internal 
motion of particles. This term is finite due to the Pauli exclusion principle. Within the hydrodynamic approach we 
have to think of this term as of an internal energy of the fluid. 

Commuting the Hamiltonian ([t]) with the density and velocity operators one obtains the continuity and the Euler 
equations of quantum hydrodynamics. Alternatively, using [fci(a;), fci(?/)] = — [A:/j(a;), fcfl(y)] = 2TriS'{x — y) the 
equations of motion can also be written as a system of quantum Ricmann-Hopf equations 



fcfl,,L + ti kRX dxknx = . (9) 
For free fermions with spin, we simply add the Hamiltonians Q written for spin up and spin down fermions: 

1 2 1 2 i-^^^ 

2PT«T + 2^i"i + -Y 



H= d.\ ^p,v^ + ^pivf + ^ (p? + pI) . (10) 



Expanding (10 1 around the background density po = — and the background velocity = up to quadratic terms 
in Va and Spa = Pa — po, 'Vfe obtain the harmonic fluid approximation 

H ^ ^ fdx {v^ + n^hHp] + vl + n^tfSpl) 
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E / [(5.0fl,a)' + (a.^L.a)'] (11) 

with right and left bosonic fields defined as dx4>ii(L),a — Va/h±TTSpa. This procedure is equivalent to the conventional 
linear bosonization procedure where the fcrmionic spectrum is linearized at the Fermi points. 
In the spin-charge basis. 

Pes = Pi±Pi and Vc^s = — ^ , (12) 



the harmonic theory (11 1 is described by a sum of two independent harmonic fluid Hamiltonians, one for charge and 



the other for spin degrees of freedom 

H^^ jdx (4i;2 + TT^nHpl + 4i;2 + n^h^Spl) . (13) 

After linearization, the quantum Ricmann-Hopf Eq. ^ reduces to (where ± stands for x — {Rj L\ respectively) 

ka^^ ± tmpa dxka,x = , a = {],[} ; (14) 

from which we identify that the quadratic excitations propagate like wave equations with sound velocities u^harge = 
Uspin = T^fipo, equal for spin and charge. Turning on interactions between fermions generally renormalizes spin and 
charge sound velocities differently and results in genuine spin-charge separation at the level of harmonic approximation. 
The spin-Calogero-Sutherland model happens to be very special in this respect. Despite a non-trivial interaction for 
spin and charge, their sound velocities remain the same. 

Although spin and charge are not truly separated for a free fermion system (and for the sCM), the interaction 



between spin and charge is absent at the level of harmonic approximation (13 1. This interaction appears if nonlinear 



corrections to ( 13 1 are taken into account (e.g., by the fully nonlinear Hamiltonian (|10[)) and due to gradient corrections 



to the hydrodynamics. The latter are not considered in this paper. 



In the proper classical limit fi. ^ all terms of ( 10 1 but the velocity terms vanish (Fermi statistics does not exist 
for classical particles). Instead, we are interested in a "semi-classical" limit in which p ^ v/h. In this limit we rescale 
time and velocity by h (t t/h and v — s- hv) and measure everything in length units. This is equivalent to dropping 



all H from equations. For instance, the Hamiltonian (10 1 becomes 



H = I dx\ ^PT«T + Ipi^I + Y (^T + pI) \ ■ (15) 
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We replace the commutation relations ^ by the corresponding classical Poisson brackets (for up and down species) 

{pa{x),vi3{y)} = SaisS' {x - y) (16) 

and consider the classical equations of motion generated by the Hamiltonian together with the Poisson brackets. In 
the remainder of the paper all hydrodynamic equations are obtained in this semi-classical limit. 



III. THE SPIN-CALOGERO MODEL 



In this work we concentrate on the hydrodynamics of sCM (jlj for the case of spin- 1/2 fermions with an anti- 
ferromagnetic sign of interaction. It is convenient to impose periodic boundary conditions, i.e. consider particles 
living on a ring of the length L. This Hamiltonian is given by 



H 



2 ^ dx"] 



A(A-P,, 

,2 TT 



sin I (xj - xi) 



(17) 



and is known to be integrabl^. All eigenstates of (17 1 can be enumerated by the distribution function 

Z.(k) =i.^(K) + i.^(K). (18) 

Here, k are integer-valued quantum numbers identifying a given state in a Bethe Ansatz description and j/|_|(k) = 0, 1 

depending on whether a given k is present in the solution of the Bethe Ansatz equations. 

The total momentum P and energy E of the eigenstate are given in terms of the distribution function iy{K) ad^fiUll 



P = 



277 



+00 



Eo + 

+00 



e= k2i/(k)-|--^|k-kXk)z/(k'), 



(19) 

(20) 
(21) 



K,— — 00 



where Eq = ^N{N'^ - 1) is the energy of a reference statd^. The numbers of particles with spin up and spin down 
are separately conserved in (17 1 and are given by 



+ 00 



K — — OO 



The ground state wave function for (17) ii 



tpGS = Yi 



jsin ^(xj - a;;)| |^sin ^(x^ - x/) 



5{aj,<7i) 



exp 



3<l 



?-sgn(crj -ai) 



(22) 



(23) 



and corresponds to the distributions 



i/|(k) = 6l(-iV|/2 < K < 7V|/2) , 
v^iti) = e{~Nj2<K<Nj2) . 



(24) 



IV. GRADIENTLESS HYDRODYNAMICS OF SPIN-CALOGERO MODEL 

Following the example of free fermions, we consider a uniform state specified by the following distributions 



(25) 
(26) 
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This state is the lowest energy state with given numbers of particles, momentum, and total spin current. It is specified 
by four integer numbers KL,fl;T,i- AH physical quantities such as energy, momentum, and higher integrals of motion of 
the state can be expressed in terms of these numbers using ( 22|19|20 1. These conserved quantities written as integrals 
over constant quantities are: 



dx Pa = — / da; 



2tt 



a = {T,i} 



"={T,i} 



da; 



47r 



Comparison with ( 2j3 1 suggests the following hydrodynamic identifications: 



, _ 277 



L 



(27) 
(28) 

(29) 
(30) 



In the main body of the paper we use vi^^i and refer to them as to "velocities" . At this point they have been introduced 
"by analogy" with the case of free fermions. In Appendices |A|B|C] we show that these velocities are indeed conjugated 
to the corresponding densities and explain their relations to the true hydrodynamic velocities. In fact, in the most 
interesting case to us, namely the CO regime (see below) these velocities coincide with the true hydrodynamic velocities 
defined in Appendix [b| The total momentum (28 1 of the system in terms of ( 29|30 ) is 



P 



da;(p|W| 



(31) 



One can also express the energy (20 1 in terms of these hydrodynamic variables. Because of the non-analyticity 



(presence of an absolute value) in formula ( 20 1 it is convenient to consider different physical regimes. These regimes are 
defined by the mutual arrangement of the supports of the distribution functions ( 25|26 1. There are six different regimes, 
that reduce to three physically non-equivalent ones using the permutation t ^ i- The distributions corresponding to 
different regimes are shown in Fig. [TJ 

• Complete Overlap (CO) regime. The support of vi is completely contained in (or vice versa). This is the 
regime considered in Ref.*, where its exact solution was given. 

• Partial Overlap (PO) regime. The supports of and of only partially overlap. 

• No Overlap (NO) regime. The supports of and of do not overlap at all. 

Notice that the small fluctuations around the singlet ground state (with ps — 0) belong to first two regimes. 

In terms of the hydrodynamic variables the three regimes are summarized in Fig. |2] and are defined by the following 
inequalities: 



Complete Overlap 
Partial Overlap 
No Overlap 



2 

TT 



<^\Ps 



Pc < \Vs\ , 



(32) 
(33) 
(34) 



where we switched to the spin and charge degrees of freedom defined by (12 1. To simplify the presentation we give 
here formulas only for the CO regime 



T^Ps 

2 



T^Ps 



(35) 



where we also assumed that p^ > 0. The opposite case < can be obtained exchanging up and down variables. 
The other regimes and formulae valid for all regimes are considered in detail in Appendix [C] 
In the CO regime ( 35 1 , the Hamiltonian can be written as 



CO 



da; 



1 



1 



Pi^l 



X 



-(/'T+/'i) + ^(2PT 



3pt/?j. 



(36) 
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(a) - Complete Overlap 



2+ 

1- 



K 



K 



K 



(b) - Partial Overlap 



2+ 
1- 



'Q.t '^Ri X'Rt 



K 



2+ 

1- 



(c) - No Overlap 






K 



K 



K 



FIG. 1: Distribution functions are shown for the three nonequivalent regimes: Complete Overlap in (a), Partial Overlap in (b) 
and No Overlap in (c). Three additional regimes exist, but are physically equivalent to the ones considered in these pictures 
and can be obtained by exchanging T ^ i- 



Complete Partial No 
Overlap Overlap Overlap 
• • • 

2 l^^l 2^'' 



FIG. 2: Diagram capturing all cases 



It is obtained by expressing ( 20|21 1 in terms of hydrodynamic variables ( 29|30 l using (35 1. As in the case of free 
fermions (see Sec. TO we now consider and V'^^i{x,t) as space and time dependent classical hydrodynamic 

fields with Poisson orackets (16 1. Of course, going from the energy of the uniform state ( 25|26 l to the nonuniform 
hydrodynamic state we neglected gradients of density and velocity fields. We refer to this approximation as to 
gradientless hydrodynamics. The equations of motion generated by the Hamiltonian (|36| with Poisson brackets ( 16 1 



can be used only when gradients can be neglected compared to the gradientless terms. This means that one can use 
this gradientless hydrodynamics only at relatively small times (compared to the time of the gradient catastrophe, see 
the discussion below). 



Before analyzing more general case let us consider some special limits of (36 1. 



Spinless limit 



In the fully polarized state = we obtain from (36) the gradientless Hamiltonian for spinless Calogero-Suthcrland 
model 



^spinless 



+00 



rfa:<! ^pw'+y (A+l)'p3 



(37) 



where we dropped the subscript f. The hydrodynamics (37) was used in to calculate the leading term of an 



asymptotics of a particular correlation function (Emptiness Formation Probability) for the Calogero-Sutherland model. 
It can be, of course, obtained by dropping gradient terms in the exact hydrodynamics derived using collective field 
theorjE™!). 
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2. A = - /ree fermions with spin 



At the particular value A = the sCM reduces to free fermions with spin and the Hamiltonian ( 36 1 becomes the 
collective Hamiltonian for free fermions (151. 



3. X —f oo limit. 



In the limit of large coupling constant X oo the particles form a rigid lattice and charge degrees of freedom 
essentially get frozen^°. We expect to arrive to an effective spin dynamics equivalent to the Haldane-Shastry modeP^'^ 
(see Appendix |d]) . This reduction to the Haldane-Shastry model is usually referred to as freezing trich^. We analyze 
this reduction in more detail in Sec. IVII 



V. EQUATIONS OF MOTION AND SEPARATION OF VARIABLES 

A. Equations of motion 



The classical gradientless hydrodynamics for sCM is given by the Hamiltonian ( 36 1 with canonic Poisson's brackets 



( 16 1 . The classical evolution equations generated by this Hamiltonian are 

Pi = -dx{pivi~ Xpiiv^ ~vi)} , 



^1 



y + ^(PT + Pi) +A7r'(P?+PTPi) + yPt 
^{vl^X. ,2^7r2A^ .9 Att^ 



(38) 



This is the system of continuity and Euler's equations for two coupled fluids (with spin up and spin down). We can 
also rewrite it in terms of spin and charge variables (12 1 



Pc 
Ps 



-dx {PcVc + PsVs} , 

-dx {Ps{Vc - 2XVs) + (2A + l)pcVs} : 



= -dx 



+ (2A + 1)^ + 



| + ^[(2A + l)V^ + (2A + l)p^]| 



(39) 



Vs = -dx { VcVs - + —ps [(2A + l)pc - \ps\ 



One can see that spin and charge are not decoupled. It turns out, however, that the variables nevertheless separate 
and the system of four coupled equations ( 38 1 can be written as four decoupled Riemann-Hopf equations (similar to 
([9])) for a special linear combinations of density and velocity fields. In the following we study the interaction of spin 
and charge governed by the above equations. 



B. Free fermions (A = 0) and Riemann-Hopf equation 



At A = equations (|38| become the hydrodynamic equations for free fermions. Fluids corresponding to up and 
down spin are completely decoupled 



n f 1 2 2 

Let us introduce the following linear combinations of densities and velocities 

kRi,Li = vi± npi. 



(40) 
(41) 

(42) 
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These combinations are nothing else but right and left Fermi momenta of free fermions. 
so-called Riemann-Hopf equation 

Ut + UUx = 0. 



All of them satisfy the 



(43) 



The equation is the same for all four combinations u — kji,L-^,i and the system ( 40|41 1 is equivalent to four decoupled 
Riemann-Hopf equations. 



The Riemann-Hopf equation (43 1 is easily solvable with the general solution given implicitly by 

u ~ uq{x — ut). 



(44) 



Here ua{x) is an initial profile of u{x, t) ai t = 0. One should solve (44 1 with respect to u and find u(.t, t) - the solution 



of (43 1 with u{x, t — Q) — uq{x). The solution (44 1 can also be written in a parametric form 



X = y + tuo{y), 
u{x,t) = uo{y). 



(45) 



This solution corresponds to the "Lagrangian picture" of fiuid dynamics and states that points in the x — u plane 
are just translated along x with velocity u, i.e., {x,uq) (x + two,uo)- This picture is especially useful to solve (43 1 
numerically. 

We notice here that the nonlinear dynamics ( |43[ ) without dispersive (higher gradient) terms is ill defined at large 
times. For any initial profile uo(a;), at large times t > tc infinite gradients Ux will develop - gradient catastrophe - and 
solutions of (44 1 will become multi- valued. The classical equation (43 1 will not have a meaning for t > tr_. We refer 



to the time tTT^unction of the initial profile) as to the gradient catastrophe time. The gradientless hydrodynamics is 
applicable only for times smaller that tc^We will discuss in more detail about validity of gradientless hydrodynamics 
in Sec. IVTIl 

We present a simple illustration of the density and velocity dynamics for free fermion system in Fig. [3] It is enough 
to consider only up-spin as the evolution of up and down spins is decoupled. We chose the initial profile of the density 
as Lorentzian with the half-width a and height h 



Pot (x) 



1 + {x/ay 



(46) 



and an initial velocity zero. We find the initial profiles of k-^-R^L using ( |42|. T hen we solve the Riemann-Hopf equations 
(43 1 using (45 1 and find the density and velocity at any time inverting ( |42[ ). We remark that for an arbitrary smooth 
bump of height h and width a the gradient catastrophe time can be estimated as tc ~ ^ . For the evolution given by 
(43 1 with an initial Lorentzian profile {uq(x) given by (46l) one can compute the gradient catastrophe time exactly. 
An infinite gradient dxU oo develops at the time 



tr 



8 a 



(47) 



For arbitrary initial conditions we compute the gradient catastrophe time numerically. 



C. Riemann-Hopf Equations for sCM 



Although the system of equations ([38| is a system of four coupled nonlinear equations, it allows for a separation of 
variables. Introducing the linear combinations of fields 



^flT.iT = W| ±7r[(A+ l)p| + ApJ = (u| ±7rp|) ± AttPc, 

kRiM = (A + l)«4-A«|±^(2A + l)p4 = («4±7rpi) + A(-2w,±2^P4) 



(48) 
(49) 



we obtain the Riemann-Hopf equation (43) separately for all four u = k^^R^-^^i. 

This property of variable separation is shared with the free fermion case Sec |VB| We notice, however, that in the 
case of sCM, variables separate only in gradientless approximation. The gradient terms neglected in this paper will 
couple the hydrodynamic equations in an essential non-separable'*'' way. 

The separation of variables in terms of ( 48|49 1 is not so surprising. One can recognize ( 48|49 ) as dressed (physical) 
"Fermi" momenta of (asymptotic) Bethe Ansatz. The integrals of motion of sCM are separated in terms of these 
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Pt(x) 



vt(x) 




FIG. 3: Dynamics of density field Pt{x) (left panel) and of velocity field v^{x) (right panel) for free a fermion case (A = 0). 
The initial density profile at f = is a Lorentzian (461 of height h — 0.25 and half- width a — 4. The initial velocity is zero. 




FIG. 4: Phase-space diagram of a hydrodynamic state characterized by four space-dependent Fermi momenta. 

Fermi momenta and the same is true for the equations of motion. We do not interrupt the presentation with this 
connection with the Bethe Ansatz solution of sCM but devote the Appendix [X] to this purpose. 

It is convenient to summarize the gradientless hydrodynamics of sCM by the picture in a "single-particle" phase 
space showing space-dependent Fermi momenta.'^ We plot the space-dependent Fermi momenta in an x — A: plane 
as four smooth lines. In the CO regime considered here (see Appendix [C| for other regimes) the Fermi momenta are 
ordered as 



(50) 



We fill the space between those lines with particles obeying the following rules of particles with fractional exclusion 
statistic^S] (see Appendix [A| (i) each particle occupies a phase space volume 2tt{X + 1) if there are no particles of 
the other species in this volume, (ii) two particles with opposite spins occupy a phase space volume 2it{ 2\ + 1) (or 
27r(A + 1/2) per particle). The velocity v^{x) is visualized as a center of a spin-up stripe on Figure |4] (see (B5 ]). The 
interpretation oivAx) is a bit less straightforward. It should be thought as a weighted average of positions of centers 
of both stripes ( |B5[ ) . 



VI. FREEZING TRICK AND HYDRODYNAMICS OF HALDANE-SHASTRY MODEL 

Here we consider the limit of large coupling constant A ^ oo. In this limit we expect that particles form a one- 
dimensional lattice and only spin dynamics is important at low energies. We refer to this limit as to a freezing of the 
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charge. We are interested in fluctuations around the uniform state with a given charge density. It can be seen from 
Fig.|4]that particles occupy the volume 2n{\ + 1/2) of the phase space when both species are present. Therefore, the 
natural expansion parameter is /i = A + 1/2 instead of A.ES We will see that the leading in term of the dynamics 
results in charge freezing, while the next to leading term gives the non-trivial spin dynamics of the lattice model 
known as the Haldane-Shastry modeP^^ 



j<i (j - 



(51) 



This model is known to be integrable.^ The freezing procedure described is referred to as "freezing trick" and was 
introduced by Polychronako^^. Our goal is to implement the procedure in a hydrodynamic description. 

Before proceeding to a regular expansion of the equations of motion we start with a heuristic argument. Wc rewrite 
the hydrodynamic Hamiltonian ( |36[ ) in terms of spin and charge variables ( 12 ) and consider first the two leading terms 
in a l//i expansion 



H 



p- 



1 J " 2 

1 Pc + p 



PcV, - PsV^ + 



2 , T^^PcPI 



PsVs 



^'pI 



9 9 

"^Pl 



-PPcPs 



12 



P 



Ps 



12 



0(1) 



(52) 
(53) 



The first term proportional to /i^ comes from the energy of a static lattice while the second term proportional to /i gives 
the Hamiltonian of the Haldane-Shastry model in the hydrodynamic formulation (see Appendix [Pj), i.e., describes the 
spin dynamics. Note that pc here should be considered as a constant equal to the inverse lattice spacing of the charge 
lattice. 



To build a systematic expansion in 1 / p we go to the hydrodynamic evolution equations given in ( 39 ) . We introduce 
the following scries in l/p = 1/(A + 1/2) for the space-time dependent fields. 



1 

— f 
P 



,(1) 



(2) 



Pc,Vc,Ps,Vs 



and re-scale time t — r/p (or dt — pdr)- We substitute these expansions into (39) and compare order by order in p. 
Let us consider few leading orders explicitly. 



A. 0{p) 



In this order the only non-trivial equation gives 



and implies that p^c^ is constant in space. 



= -9. 



(54) 



At this order we have 







B. 0(1) 






= 0, 










;2p(").(")-2p(°)«f 




,(0) 








,(0) 









(55) 
(56) 

(57) 
(58) 



spin density, ps^^ and ( 58 1 for the spin velocity Vs"' do not depend on the dynamics of the charge and are precisely 



Combining (54 1 and (55 1 we see that pi*^^ is a constant independent of space-time. The evolution equations (56 1 for 



,-,(0) 
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the ones obtained for the Haldane-Shastry model (compare to (D17l). We refer the reader to the Appendix [P] for 
more details on the hydrodynamics of the Haldane-Shastry model. 



Equation ( 57 1 is important in resolving a well-known "paradox" . In the original spin-Calogero model the momentum 



of the system is identical to the total charge current since all particles in the model have the same charge. On the 
other hand in the Haldane-Shastry model the momentum is carried by spin excitations and superficially no charge 
motion is involved. One can ask how this is compatible with getting the Haldane-Shastry model in the limit A ^ oo 
from the spin Calogero model. Equation (57 1 is necessary to make sure that the the current density jjx) = PcVcA- Ps'^s 
is globally conserved at a give n order in 1/ p,. Since pi'^ is a constant in space-time we expect from (56) and (58) that 
v'^^ evolves according to (57) to ensure that the current density is conserved. As a result, there is a charge motion 



associated with the momentum but in the large A limit this "recoil" momentum is absorbed by the whole charge 
lattice. 



C. 0(1/m) 



For the sake of brevity we do not write down the equations at this order but make some comments instead. In 
the previous order, 0(1) we noticed (see eqs. (56) and (58) that spin degrees of freedom evolve as the charge is 
essentially frozen and at that order there is no feedback of the charge degrees of freedom on the spin. However, 
in the order 0(l//i) we have feedback terms in both evolution equations for and Vs- As an example we have 

pi^^ = -9. {... + 2vf\^P + vi'^pf^ + ...} and vi'^ = -9. {... + v^'K?^ + ^ pi^Vi"^ + ••■ 
there is a charge feedback into the spin sector. 



. I which clearly show that 



D. Evolution equations for Haldane-Shastry model from the freezing trick 



The shortest way to evolution equations for Haldane-Shastry model is to take A oo limit directly in Ricmann-Hopf 
equations ( 43 1 . After rescaling time t — r/p, we have 



where k = kj p = fc/(A -I- 1/2). In the large A limit we have using ( 48|49 1 fc_R|,L| — > i'l'Pc and kuxxi 
Then the equation (59) gives evolution equations for Haldane-Shastry model with ( D12|D15 ). 



(59) 



VII. ILLUSTRATIONS 



It is relatively simple to obtain the evolution of arbitrary (smooth) initial density and velocity profiles solving 
equations of the gradientless hydrodynamics ( 39 1 numerically. One can do it very efi^ectively using the fact that the 



dynamics is separated into four Riemann-Hopf equations (43) and using their general solutions (45). In this section 



we give numerical results for charge and spin dynamics corresponding to a relaxation of a (spin) polarized center. 
These results show that due to nonlinearity of equations spin can drag charge in spin-Calogero model. We notice here 
that in the examples considered in this section the dynamics belongs to CO regime.'^ 



A. Charge dynamics in a spin-singlet sector 



As a first example we consider initial conditions ps-,Va = ^ and some arbitrary initial conditions for pc and v^. It 
is easy to see from ( 39 1 that spin density and spin velocity remain zero at any time while charge degrees of freedom 
satisfy 



Pc = ~dx{PcVc), 




(60) 



Hydrodynamics ( |60[ ) is identical to the one of the Calogero-Suthcrland model with one species of particles (except 
for the change A+l^A+1/2). It can be written as a system of two Riemann-Hopf equations (43) for fields 
± 7r(A + l/2)pc. 
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We conclude that the charge dynamics does not affect spin in a spin-singlet state at least in the gradientless limit. 
It is interesting to see how spin dynamics affects the charge one. 



B. Dynamics of a polarized center 

To see how spin drags charge we start with an initial configuration with static and uniform charge background. We 
assume that initially there is no spin current but there is a non-zero polarization given by a Lorentzian profile: 

h 

i.e., there is an excess of particles with spin up over particles with spin down near the origin. The maximal polarization 
is h and a half-width of the polarized center is a. As an illustration of spin and charge dynamics we present a solution 



of (39 1 with initial conditions (61 1 corresponding to h — 0.25 and a ~ 4. Some important comments are in order. 



Applicability of gradientless hydrodynamics 



The hydrodynamic equations we use (39) neglect gradient corrections and, therefore, are approximate. They can 



be applied only under the condition that the neglected higher gradient terms are small compared to the terms taken 
into account in (39). Of course, the exact criteria can be written only when the form of the higher gradient terms 
are known explicitly. Here, we are going to use a much simpler criterion. We require that all fields change slowly 
at the scale of the inter-particle spacing. The uniform background pc — I defines the inter-particle spacing and the 
characteristic scale for hydrodynamic fields to be 1 and we require d^f I for all fields at all x and t that we consider. 

One can easily check that dxPs{x, t = 0) <C 0.1 for all x with the initial profile (61 1 (in fact, the maximal derivative 
is approximately 0.041). Because of the gradient catastrophe this condition will be broken at some time and we 
can trust the results obtained from (39) only up to that time. To be well within this criterion all our fields satisfy 
dxf < 0-3 at any given time. 

Let us start with the solutions for the case of free fermions, i.e., A = 0. 



2. Free fermions with spin: A = 



We present the results for spin and charge dynamics of free fermions with polarized center initial conditions (61) 
on left panels of Figures 5]6 The profiles Ps{x) and pdx) — 1 are shown as functions of x for times r = 0, 1, 3.5, 7 



respectively. Here we use a rescaled time r = (A -I- l/2)t = t/2 for future convenience. 

The dynamics is separated into four Riemann-Hopf equations for each Fermi momenta. The initial conditions (61 1 
can be written as Lorentzian peaks for each of the four Fermi momenta of fermions and all four Fermi velocities are 
different. This results in a splitting of an initial Lorentzian peak into four peaks at larger times which can be easily 
seen on the left panel of Fig. [5j In addition to this linear effects the nonlinear effects of steepening the wave front can 
also be seen. The latter will render gradientless hydrodynamics inapplicable at later times. 

The drag of charge by spin clearly seen in Fig. |6] has an essentially nonlinear nature. There is an excess (deficit) of 
particles with spin up (down) at the origin at the initial moment. The particles with spin up will move away from the 
center while spin down particles will move towards the center. However, the average velocity of spin up particles is 
larger than the average velocity of spin down particles as it is proportional to the density of those particles. Therefore, 
the initial motion of particles away and towards the origin creates a charge depletion in the center and charge density 
maxima away from that depletion. This gives a qualitative explanation of the picture of charge dragged by spin which 
is shown in the left panel of Fig. |6] Notice that in this explanation we used the dependence of propagation velocity 
on the amplitude of the wave - an essentially nonlinear effect. 



3. X-dependence of spin and charge dynamics 



To see the effects of the interaction on spin and charge dynamics we show the spin and charge density profiles 
at a fixed time for different values of the coupling constant A in the right panels of Figures |5|6| respectively. It is 
convenient to use the scaling dictated by the A — *■ c» limit considered in detail in Section |VI[ Namely, we use a 
rescaled time r = (A + 1 /2)< and rescale the deviation of the charge density from the uniform by background plotting 
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FIG. 6: Left panel: Charge dynamics of polarized center for free fermions. The initial charge density profile is constant and the 
initial spin density profile is a Lorentzian ( 61 1 of a height h = 0.25 and a half-width a = 4. Profiles at times t = t/2 — 0,1, 3.5, 7 
are shown. Right panel: A snapshot of a rescaled charge density (A+l/2)(pc — 1) at time t = r/(A-|-l/2) for t — 7 for A = 0, 1, cxa. 



(A + 1/2) (pc — 1) for the charge density. The charge and density profiles found at r = 7 are remarkably close for A 
ranging from the free fermion case A = to the limit of Haldane-Shastry model A — > oo. 

The results confirm that the effect of spin dynamics on charge is suppressed by 1/A for large A. For a given initial 
spin density profile the maximal amplitude of charge deviation is of the order l/(A + l/2). 



VIII. CONCLUSIONS 



In this paper we considered a classical two-fluid hydrodynamics derived as a semiclassical limit of the quantum 
spin-Calogero model (sCM) defined in (17 1. The model (17 1 is essentially quantum as it involves identical particles 
and a particle permutation operator. There is an essential ambiguity in how one takes a "semiclassical" limit. Here 
we considered a limit which is obtained when the density of particles goes to infinity so that hp is kept finite in the 
limit h 0. We have also neglected gradient corrections to hydrodynamic equations assuming that fields change very 
slowly on the scale of the inter-particle spacing. With all these assumptions, hydrodynamic equations are obtained 
from the Bethe ansatz solution of sCM. They have the simplest form when written in terms of fields corresponding 
to dressed Fermi momenta of Bethe ansatz. In terms of these fields ( 48|49 ), equations separate into four independent 
Riemann-Hopf equations (43 1 which are trivially integrable. 

We presented some particular solutions of the hydrodynamic equations illustrating interactions between spin and 
charge. There is no true spin-charge separation in sCM. However, in the limit of large coupling constant A ^ oo the 
spin degrees of freedom do not affect the dynamics of charge degrees of freedom. The spin dynamics then is described 
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by the hydrodynamics of Haldane-Shastry spin model. We considered expHcitly both this hmit (A — > oo) and the 
hmit (A = 0) of free fermions with spin. 

The quantum scattering phase of particles interacting via potential is momentum independent. Moreover, 
it is the same for particles of the same species and for particles of different species because of the SU{2) invariance 
of (17 1. It is well known that this allows one to describe sCM as a model of free exclusons - particles obeying an 
exclusion statisticEpmnilol. We do not keep the SU{2) invariance of the original quantum model il7h explicitly when 
taking the classical limit. However, this invariance is responsible for the variable separation that we observed in our 
hydrodynamics. We note here that sCM can be generalized to the "multi-species Calogero model"^. Because of the 
absence of the SU{2) invariance for a more general two-species Calogero model one does not have the separation of 
variables for the corresponding hydrodynamics. 

The classical gradientless hydrodynamics derived in this paper captures a lot of the features of sCM. It is straight- 
forward to generalize our results to the case of the SU{n) Calogero model and to use the gradientless hydrodynamic 
equations for problems where field gradients can be neglected. In a separate public at ion^SI we use these equati ons in 
instanton calculations for the computation of emptiness formation probability similar to what was done in Ref^^MH 

However, some important features of the hydrodynamic description do require an account of gradient corrections. 
First of all, the exact hydrodynamic equations are expected not to have an exact separation of variables. The 
obtained Riemann-Hopf equations ( |43[ ) acquire gradient corrections and four such equations written for ( 48|49 1 are 
expected to be coupled by those gradient corrections similarly to the case of the one-species Calogero-Sutherland 
modepi. Similarly, we expect that the equations w ith gradien t corrections will have soliton solutions corresponding 
to quasi-particle excitations of the quantum model ( I7) 6 | i5 | 34 ]^ 

The hydrodynamic description of quantum sCM has been addressed in Refs.'^^^ using the collective field theory 
approach. The comparison of our results with results of those works is not straightforward. One should apply the 
collective formulation of Refs.'^^'i^ to the states from an appropriate sector of coherent states and take a corresponding 
classical li mit. It would be especially interesting to see how the three hydrodynamic regimes discussed here appear 
from Refs.li^. One can also recognize a lot of similar looking terms in quantum hydrodynamics of Refs.'^ and in 
our classical gradientless hydrodynamics. It would also be very important to understand the role of the degeneracy 
due to the Yangian symmetry in sCM on its hydrodynamics. The latter degeneracy was neglected in the classical 
hydrodynamics of this paper. 
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APPENDIX A: ASYMPTOTIC BETHE ANSATZ SOLUTION OF SPIN CALOGERO MODEL AND 
SEPARATION OF VARIABLES IN HYDRODYNAMICS 



The spin Calogero model is solvable by asymptotic Bethe Ansatz (ABA^^I^MI] This solution turns out to be the 
most convenient for our purposes. 

The most important ingredient of ABA is the scattering phase which is given by 



e{k) = 7rAsgn(fc) 



(Al) 



for sCM. Here k is the relative momentum of two particles and the scattering phase does not depend on the species 
of particles. The expression for the dressed (true physical) momentum of the particle is given by 



fc(K) 



2tt 

T 



sgn(K — k')v{k') dn! 



(A2) 



where k is an integer- valued non-interacting momentum of the particle (quantum number) and v{k) is the number of 
particles with quantum number k (see (18)). Here we replaced in the scattering phase sgn (fc — k') by sgn (k — n'^^. 



16 



We immediately obtain from (A2| 



dfc 

L— = [1 + \iy{K)] dn 

ZTT 



and 



^T(i) 



(k) dn = L 



V](i){k) dfc 



1 + \v{k) 27r' 



(A3) 



(A4) 



We can see that the picture corresponding to ( |A4| in a single-particle phase space requires that the number of particles 
in the phase space volume is given by 2^(\+i) °^^y "-"^"^ species is present and 2tt(\+i/2) ^^^<5n both species are present. 
This justifies the picture we used (see Figs. |4]and[7]). 

It is easy to write down the expressions for the conserved quantities using ( A4 1 . 

^+^dk 



T(i) 
P 

Ps 

E 



— OO 
+ 00 



— OO 
+ 00 



— OO 
+ 00 



2tt 


1 


+ \u{k) 


dfc 




v{k) 


2tt 


1 


+ \y{k) 


dfc 




Vs{k) 


2tt 


1 


+ \y{k) 


dfc 




v{k) 


2^ 


1 


+ Xv{k) 



Here Pg is a conserved quantity proportional to C\ introduced i: 

N 

2 L 



(A5) 
(A6) 
(A7) 
(A8) 

(A9) 



One can think of, e.g., P| = (P + Ps)/'^ as of a sum of asymptotic values of momenta of spin-up particles. We have 
replaced summations by integrations as we need only continuous versions of these formulae. It can be shown that 
( A5 A6|A8 1 are equivalent to ( 22|19|20 ) with the relation between physical and non-interacting momenta given by 



(|A2|). Moreover, because the measure of integration |^ i+\j{k) ^® * piece- wise constant for the two-step distribution 



one naturally obtains integrals of motion in a form which is completely separated in terms of Fermi momenta. 
Indeed for a two-step distribution ( 24 1 



1, if k^a < k < k 
0, otherwise 



Ra 



(AlO) 



where a =t, i- In the CO regime (50 1 we have 
dfc J^t(fc) f'^'-i dk 



2tt l + Xi^ik) 

^ '^i(fc) 
27r 1 + \v{k) 



f{k) = 

m - 



1 



27r IH- A 
dfc 1 



27r 1 + 2A 



fik) + 
f(.k), 



dfc 1 



27r l-f 2A 



fik) 



dfc 1 



27r 1 + A 



f{k), 



where /(fc) is an arbitrary function. In particular, we obtain for the densities 

kL^ + 2A + 1 *^^^^ ~ 



27r(A- 






27r(A + 




= kn^ 


47r(A- 


-•f 


- fc2 


47r(A-f 




- fc2 


127r(A- 




- fc3 



2A + 1' 
1 



fc3 



2A-h 1 



kii)- 



(All) 

(A12) 
(A13) 
(A14) 
(A15) 
(A16) 
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So far we presented the values of the conserved quantities for the sCM in terms of dressed Fermi momenta. They 
are given by hnear combinations of Fermi momenta raised to the same power. There are infinitely many integrals 
of motion of this type and they are all in involution (commute with each other). The latter is a pretty stringent 
requirement and we assume that the only way to satisfy it is to require that the corresponding classical hydrodynamic 
fields have the following Poisson's brackets 



{ka{x),kp{y)} = 2TrSaSaf3S'{x - y), 



(A17) 



where a runs over all Fermi points and Sa are some numbers to be determined. We can determine these numbers, 
e.g., in the following way. The density of current j (momentum per unit length) from (A 14 1 by 



1 



47r(A + 1) 



"■AT 



1 



2A + 1 



(A18) 



The total momentum of the system is a generator of the translation algebra {P,q(y)} — dyq{y), where q{y) is any 
field. For the current density we should have 



{ji.x),q{y)} = q{:x)5'{;x - y). 



(A19) 



Taking q(y) to be ka{y) and combining (A19l with (A17) we fix the unknown coefficients s 



Sri 



-sl-\ = a + 1, 
-SLX-(A + 1)(2A + 1). 



Computing Poisson's bracket of the hydrodynamic Hamiltonian (obtained from (A16l 

1 



H 



127r(A + 1) 



da; 



1 



2A+ 1 



with ka{x) we obtain Riemann-Hopf equation (43 1 for every Fermi momentum field ka{x,t) 



(A20) 



(A21) 



APPENDIX B: HYDRODYNAMIC VELOCITIES 



In appendix |A] we did not use the notion of hydrodynamic velocity. Instead, our hydrodynamic equations were 
written directly in terms of dressed Fermi momentum fields ka{x,t). We also know how to express other quantities 
like density, momentum, energy, etc in terms of these variables. Let us now find the expressions for the velocity fields 
We focus on the CO regime here and consider other regimes in appendix [C| 

First of all we, give the expressions for the conserved densities and conserveocurrent densities which can be found 
from ( |A12|A13|A14|A15| ) as 



Pi 

it 
h 



p 


+ Ps 


1 






2 


27r(A4 


1) 


p 


- Ps 


1 






2 


27r(2A 


+ 1) 


3 


+ 3s 


1 






2 


47r(A + 


1) i 


3 


-3s 


1 






2 


27r(2A - 


fl) 



kR] 


- kL] - 


-{kRi 


-kLi). 


k'^ 


k^ 


(kill 


-klO- 



A 



2A 



^{kRi - kLi) 



2A+ 1 



(4i 



klO 



(Bl) 



(B2) 



In hydrodynamics, the velocities are defined as variables conjugated to the conserved momenta. Namely, the differ- 
ential of the energy density defines chemical potentials and velocities as 



de = /i|dp| + fi^dpi + v'^dj^ + v'ldj^. 



Using the energy density obtained from ( A21 1 we have 

1 



de = 



47r(A+ 1) 



2A 



- (fclj, dkRi - kl^ dkLi) 



(B3) 



(B4) 
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and using (|B1|B2D one can determine and The hydrodynamic velocities are given by linear combinations of 
Fermi momenta - 

< = 2(aVi) [^(fcfiT + fc^T) + (fcfli + fcLi)]. (B5) 

Using ( |A17|A20| one can check that the velocities te5k have canonical Poisson's brackets with densities (^BlP^I 

{pa{x),vp{y)} = 6ai3S'{x ~y), (B6) 

where a,/3 ^tii- The other Poisson's brackets vanish. 

The hydrodynamic velocit ies (B5l are precis ely the ones used in the main body of this paper for CO regime 
I = W| |. Equations ( 48|49 1 are the inverse to (Bl B5l. Interestingly, in the CO regime the velocities and densities 

of different species can be naturally (simply) written in terms of bare non-interacting momenta ( 29|30 l. 
The current density in terms of densities and velocities follows from (Al8l (compare with (311) 



(B7) 



The density of "spin-current" which follows from (A15) has a "correction" proportional to A compared to the case of 
free fermions 



(B8) 



In this appendix we focused on CO regime. Of course, the formalism reviewed here is applicable to all three 
hydrodynamic regimes (CO, PO, and NO). We collect appropriate results in Appendix [Cj 

APPENDIX C: HYDRODYNAMIC REGIMES FOR SPIN-CALOGERO MODEL 



Depending on the relative order of four quantum numbers kj^^l-i^^i we distinguish six different hydrodynamic regimes 
of the sCM. These regimes can be reduced to three essentially different ones exchanging t ^ i- Iii this appendix we 
consider these three regimes and then combine all six cases. 



Before we proceed, let us remark that the function k^n) defined in ( A2 1 is monotonic and the order of the quantum 



numbers kr.l:].i is the same as the one of the physical dressed momenta fc/j,| — fc(K_R.|) etc. Therefore, we can use 
the latter to define hydrodynamic regimes instead of the bare momenta k. 



Conserved densities and dressed Fermi momenta 



Let us consider generally some integrable system which has two infinite families of mutually commuting conserved 
quantities. We assume further that the densities of these quantities are given in terms of four dressed Fermi momenta 
ka{x) with a = 1, 2, 3, 4 as 



1 4 

jn{x) = - aa(fca(a;))", 
n ^-^ 

fn^^) = - Vac.&a(fca(a;))". 



(CI) 



a=l 



Here n — 1,2,3,... and constant coefficients. We assumed that the conserved densities can be expressed 

locally in terms of and neglected gradient corrections. 

We identify the first several integrals with densities, currents, and the energy as 



h{x 
h{x 
jsix 



= Ps{x), 
= js{x), 

= 2e{x). 



(C2) 



19 



TABLE I: Summary of three regimes. 





a 


LT 


R T 


LI 


Ri 




CO 




1 

— 1 

1 


1 
i 

1 


1 

2A+1 

-(2A+ 1) 


1 

2A + 1 

2A + 1 


fciT < ^Li < kni < kn-j 


PO 


2n{X + 1) ttQ 


-1 
1 


1 

2A+1 

2A+ 1 


1 

2A+1 

-(2A+ 1) 


1 

-1 


kLi < kLT; < kni < kn^ 


NO 


2n{X + 1) aa 

ba 


-1 
1 


1 
1 


-1 
-1 


1 

-1 


khi < kni < kLt < kn-i 



We notice here that due to ( A5|A8 ) t he i dentifications (C2| (with (Cll) are vahd for sCM model in all its regimes. 
The higher order conserved densities (CI) correspond to conserved quantities of sCM introduced in Ref."^. 



The requirement of vanishing Poisson's brackets between conserved quantities is very restrictive. It can be resolved 
by requiring canonical Poisson's brackets between Fermi momenta (A17l. If (A17l is valid, it is easy to ch eck that 
{/ dxjn{x),J dy j^{y)} = etc. Using the fact that the total current is the generator of translations (Al9l we can 
fix the coefficients Sa in (A17l as 2TrSa = ^/cia and obtain 



1 



{ka{x),kf3{y)} = — SapS'ix - y) 



(C3) 



Using the Poisson's brackets (C3l and the Hamiltonian H = J dxe{x) with ( C2|C1 ) it is easy to obtain the Riemann- 
Hopf evolution equations for the dressed Fermi momenta 



dtka + k^dxka — 0, for a = 1, 2, 3, 4 



and the evolution equations for all conserved densities as 



(C4) 



(C5) 



In the hydrodynamic regime only four of the densities are algebraically independent (as there are only four dressed 
Fermi momenta). Therefore, one can find constitutive relations, i.e., express the energy density in terms of p,ps,j 
and js. Alternatively, one can use hydrodynamic velocities and defined by (B3l instead of currents j,js- 

We can see that the hydrodynamics ( Cl|C2|C3 l is fully defined by coefficients aa,ba- In fact, these coefficients are 
not totally independent. Requiring that densities p and ps have vanishing Poisson's brackets with themselves and 
with each other gives three relations between the coefficients 



^ baao 



0, 
0, 



= 0. 



(C6) 



For CO, PO, and NO regimes of sCM these coefficients are summarized in the Table |T] These coefficients do satisfy 
relations (C6l. 



The matrix of Poisson's brackets of the dressed Fermi momenta ka (C3l is diagonal but not proportional to the 



unit matrix. It is interesting that the Poisson's brackets of bare momenta satisfy 



{Ka{x),Kfl{y)} 



{-ir^S^pS'{x-y). 
Zn 



(C7) 



One then obtains that the "velocities" introduced in ( 29|30 ) are canonically con jugate to the corresponding densities 
and can be written as linear combinations of Ka (and of ka). The velocities ( 29|30 l are defined just as conjugate 
variables to the densities. This definition is not unique. One can always shift W| — > W| + 2'K^pi and vi ~> f x — 27r7P| 
with any number 7 without changing Poisson's brackets. The particular choice of variables ( 29|30 l is convenient 
because it defines velocities continuously across all hydrodynamic regimes. Moreover, we have 



-'T.i 



± TlXp 



for CO, 
for NO. 



(C8) 
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In PO regime the hydrodynamic velocities are not linear combinations of ka and their relations to the conjugated 
variables used in this paper are more complicated. 



2. Complete Overlap Regime (CO) 



The Complete Overlap regime corresponds to the case when 

TT TT 



(C9) 



In this case the support of vi is a subset of the support of vi (or vice versa) . In the main body of the paper we mostly 
concentrated on this case, but for convenience we recap the main formulae in this appendix as well. The dressed 
momenta (A2) in the CO regime for ps > 0, i.e., for the ordering 



i?.T : 



(CIO) 



are 



fci?.T,iT = ft ± TT [(A + l)p| + Apx] = U| ± 7rp| ± AttPc 



(Cll) 



Poisson's brackets of ka are given by (C3l with coefficients from the Table |T] One can express all conserved densities 
(Cll in terms of dressed Fermi momenta using the Table |l] For example, the Hamiltonian (see (C2|) reads 



1 



CO 



127r(A+ 1) 
' 1 



dx 



k^ 



2ATT(''^^ 



(C12) 



= / dx 



^2^2 ^2 ^^2 



(C13) 



The evolution equations are given by ( C4 1 and can also be recast in terms of equations for densities and velocities 
p8p9|. 



3. Partial Overlap Regime (PO) 

There are two regimes when the supports of and v-^ only partially overlap. Here we concentrate on the case for 
which 

I Ip.I < ^'s < I Pc , (C14) 

corresponding to the ordering 

kLi < kL] < fcfli < fcfl,| . (C15) 

The other PO regime can be obtained by exchanging up and down particles, i.e. by changing Vs — > —Vs- In this case 
the dressed momenta ( A2 1 are 

kLi = - 7r(A + l)pi - TrApi = - npi - Xirpc, 

fc^i ~ vij + A(ti| — vi) ~ tt{2X + l)p| = W| — 7rp| + A(2?;s — 27r/9|), 

kRl = - - + ^^(2^ + = «i + T^Pl - H'^Vs - 27rp|), 

fcflT = ''^T + ^('^ + -'^)'^T + "^-^Pi = + 7rp| + Xnpc (C16) 



21 




FIG. 7: Phase-space diagrams of a hydrodynamic states characterized by four space-dependent Fermi momenta in three regimes 
CO, PO, and NO respectively. 



and the Hamiltonian becomes (see Table [l] and ( C1|C2 )) 



PO 



127r(A+ 1) 



da; 



^ii ~ ^LT + 2A+ 1 ^^^^ ~ ^^^-^ 



(C17) 



TT A 



6 



6 



(C18) 



Poisson's brackets of are given by (C3l with coefficients from the Table |T] and evolution equations are given by 
(p|. 



4. No Overlap Regime (NO) 

In this case, the supports of and do not overlap at all. For > the ordering of dressed Fermi momenta is 



and momenta themselves are 



kRi,Li = - TrApi ± 7r(A + l)p| = wj. ± Trp^ - Xnps^c- 



and the Hamiltonian becomes (see Table |T] and ( |C1|C2 1 ) 



127r(A-t- 1) 

dxj ^piwj -f Attpip; (w| - wi) 

2 2 

+ y (A + 1)' (PT + Pi)' - Y (1 + 2A) PTPi (PT + Pi) 



(C19) 



(C20) 



(C21) 



(C22) 



Pois son's brackets of ka are given by (C3l with coefficients from the Table |l] and evolution equations are given by 



22 



k inequality 




Ps 


6 = + f 


6 = «s - f Ps 


Xl = + f Pc 


X2 = ys- f Pc 


Vs inequality 


Regime 


fcit < fcflT < f^Li < kni 






- 


- 


- 


- 


Vs < -^pc 


NO 


feiT < f^Li < fcflt < kRi 


- 




- 


- 


+ 


- 


-^Pc <Vs < -f IPs! 


PO 


khi < fcii < fciji < kni 








+ 


+ 




f ps < Vs < —\ps 


CO 


kLf < kLi < kni < kRf 




+ 


+ 




+ 




— \ps < Vs < f Ps 


CO 


kn < kL] < kni < fcfli 


+ 




+ 


+ 


+ 




f jpsj <Vs < Ipa 


PO 


khi < kni < kL] < fc/it 


+ 




+ 


+ 


+ 


+ 


^Pc < Vs 


NO 



TABLE II: Classification of different regimes: + indicates that the field takes positive values, — that it is negative. A blank 
means that its sign is arbitrary. 



5. All cases combined 



It is possible to combine all hydrodynamic regimes into relatively compact expressions introducing absolute values 
of hydrodynamic fields. A general Hamiltonian valid for all regimes takes a form 



H 



dxl ^PT^? + Ipi^! + y (^T + p1) + y ^'/'c + y2A(p? + 



+A p.Ci6 - ^ (161' + m + ^ [Ixil' - X? + \X2\' + xl] 



where we introduced the following notations 



6,2 
Xl,2 



Vs±^ Pc ■ 



(C23) 

(C24) 
(C25) 



The Hamiltonian (C23) can be obtained from ( 20|21 ) for the general case of a two-step distrib ution function 

We collect in the Table [ll| the information necessary to go quickly from the general expression (C23 1 to the particular 

ones valid in separate regimes (CO, PO or NO). 



We can combine the evolution equations following from (C23l in the spin/charge basis (12 1 as 

Pc = -dx^PcVc + PsVsY 

Ps = -9a;< PcWs + PsUc - - kilCil + 6161 - XilXil - X2IX2I 



(C26) 
(C27) 



(4A2 + 2A + 1)p^ + (2A + 1)p^ 



XilXil - X2IX2I 



(C28) 



dxi VcVs + — (2A + 1) pcPs 



Ar 
2 



6161 -6161 



(C29) 



For CO and PO regimes the Hamiltonian (C23) takes an especially simple form in terms of dressed momenta 

A 



Ho 



1 



o&po — 



127r (2A + 1) 



(A + 1) 



Fit '^Ll\ 



which are related to density and velocity fields as 

fcflT,iT = «T =t ^(A + 1)Pt + '^Xi,2 T A|6,2| , 
I^Rl,Ll = «i ±7r(l + A)px - Ax2,i T A|6,2| ■ 



(C30) 



(C31) 



As in the separate cases considered before, these momenta have canonical Poisson's brackets (A17) with 

SRTMl = (A + l)[A + l±Asgn(6)] 

SLT.Li = -(A + l)[A+lTAsgn(6)] (C32) 



and evolve independently according to the Riemann-Hopf equations ( C4 1 . 
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APPENDIX D: HYDRODYNAMIC DESCRIPTION OF HALDANE-SHASTRY MODEL FROM ITS 

BETHE ANSATZ SOLUTION 



The Haldane-Shastry model (HSM) is a Heisenberg spin chain with long-ranged interaction defined by the Haniil- 
tonian: 



■"HSM = 9 2^ 
3<l 



where we Kji is the spin-exchange operatoil^ 



(Dl) 



(D2) 



and d{j) = (N/tt) |sin (7 rj/jV )| is the chord distance between two points on a lattice with N sites and periodic boundary 
conditions. The model (Dll has been introduced independently at the same time by Haldanft^fiO and by ShastrjED and 



has been shown to be integrable. The energy spectrum of the HSM is equivalent to that of the Calogero-Sutherland 
model at A = 2, but with a high degeneracy due to the Yangian symmetrj^^siioj^ 

In this appendix we used the Bethe Ansatz solution'^'^ ■^^ ^ to construct a gradientless hydrodynamic description for 
the HSM similarly to what we have done for the sCM model in section IV and appendix]^ To this end, we consider 
a state with M overturned spins over an initial ferromagnetic configuration (say from up to down and M < N/2) and 
introduce M integer quantum numbers k's to characterize the state in the Bethe Ansatz formalism. As before such 
state can be described by a distribution function v^k) = 0, 1, depending on whether that quantum number is present 
or not in the BA solution. FoUowing'^^ we impose a condition on the integer numbers: |k| < {N — M — 1)/2P^ 

The scattering phase for the HSM is 



which corresponds to setting A 



0{k) = TT sgn(fc) , (D3) 

1 into (AlP^l, Please note that since we are considering a lattice model, the 
-TT < fc < TT, where we took the lattice spacing as unity. 



momentum is defined within the Brillouin zone: 

At this point, all the derivations of appendix|^can be repeated step by step for the HSM just by setting everywhere 
A = 1, and remembering that the momentum is always defined modulo 27r. In particular, the dressed momentum is 



k{K) = 



2n 



K+- / sgn (k — K')l>(K')dK' 



(D4) 



where again we replaced sgn (fc — fc') by sgn (k — n'^j^ and the distribution of the physical momenta is given by 



N 

i'{K)dK — — i/(fc)dfc . 
In terms of this distribution function, the conserved quantities can be written as 

/dfc 

P ^ N f — v{k) fc, 
J 47r 

where the momentum is defined only modulo 27r. From now on, we will drop the constant energy shift Eq. 
In a hydrodynamic description we assume a distribution of the uniform type 



(D5) 

(D6) 
(D7) 
(D8) 



.(fc) = 



1, if— 7r<fcL<fc<fc_R<7r, 
0, otherwise, 



(D9) 



where fc^j^L are some numbers. Using (D9l and introducing space dependent fields instead of constants we write 
( |D6|D7| ) as 



M 



ax = ax p, 



i-2 i.2 



— I dxpv, 



(DIO) 
(Dll) 
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which suggests the identification 



Then the hydrodynamic Hamiltonian follows from (D8l 



^3 _ ^3 

dx = / da; 

247r 



(D12) 



(D13) 



which corresponds, as expected, to the (gradientless) hydrodynamic of a A = 1 spin-less Calogero-Sutherland model 
([371. 



We think of slowly varying fields p{x, t) and v{x, t) as of classical fields obeying the Poisson relation {p{x),v{y)} 
S'{x — y). Then (D13) generates the evolutions equations 



P 

i) 



-dx (pv) , 



(D14) 



One can easily recognize in (D14l the hydrodynamics of spinless Calogero-Sutherland model (37) for A = 1. The corre 



spondence between eigenstates and eigenenergies of Haldane-Shastry model with A = 2 spinless Calogero-Sutherland 
model has been noticed in the original papeif^. The degeneracy of the states due to the SU{2) invariance and Yangian 
symmetry is lost in our classical hydrodynamics model. 

For comparisons with the derivations from freezing tricliPSl in section 



VI 



we express ( |D13|D17 | in terms of ps and 



Vb used in the main body of the paper. We identify the density p — M /N — pi as the density of spin-down particles 
and the velocity w as a velocity of spin-down particles relative to the static background of spin-up particles, i.e., 
V = vij — vi = —2vs. The charge density corresponding to the lattice with spacing one is just pQ = I. We summarize 



Pi 



Po 



-2v, 



Po 



1. 



(D15) 



Using (D15) we rewrite (D13l as 



HSM 



^ i dx{ pov^ 



PsVs 



2 2 



12 



(D16) 



where we neglected a constant and a term linear in p^, which amounts to a shift in the chemical potential. The 
evolution equations for the spin density and spin velocity follow from ( D14|D15 1 



= -dx {2vsPo - 2vsPs} , 
^2 



-poPs 



(D17) 



— 2 



and 



We notice that the above (D16 D17) is nothing but the strong interaction limit of the sCM ( 53|56|58" l 
Finally, we remark that it is easy to check that the distribution function (D9) implies that — < 
therefore corresponds to the CO regime of spin-Calogero model. 

Both in this appendix and in writing classical hydrod ynami cs for sCM we neglected the degeneracy of the corre- 
sponding quantum models due to the Yangian symmetrj^^Hin!^ \Ye assumed that during the evolution string states 
are not excited. Of course, the degeneracy plays a very important role for perturbed integrable systems and for the 
hydrodynamics at finite temperatures. 
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We notice here that for a free fermion system it is possible to assign the meaning even to the multi- valued solution of 1 44 1 



for t > tc. It is the boundary of the support of the Wigner distribution in the one-particle phase space. In this paper we 



restrict ourselves to times less than the time of gradient catastrophe and assume that ( 43 1 has a well-defined single- valued 
solution. 

At least one will not be able to separate variables considering simple linear combinations of fields. To the best of our 
knowledge variables in sCM do not separate or, at least, an appropriate change of variables has not been found yet. 
We would like to thank A. Polychronakos who encouraged us to present this picture. 
Of course, first orders of the expansion are not sensitive to this shift. 

In exotic cases involving boundaries between CO and PO regimes one notices singularities developing at the boundary and 
we expect gradient corrections to correct these singularities. 
Note that for fermions P^i Kji = —1. 

This corresponds to having a single compact support of u within a single Brillouin zone. Other regimes will require an 
analysis of umklapp processeP^^ and will not be considered here. 

This scattering phase is identical to the one in A = 2 bosonic Calogero-Sutherland modeP^. 



